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We discuss the effect of electron-electron interactions on the static polarization properties of 
graphene beyond RPA. Divergent self-energy corrections are naturally absorbed into the renor- 
malized coupling constant a. We find that the lowest order vertex correction, which is the first 
non-trivial correlation contribution, is finite, and about 30% of the RPA result at strong coupling 
a ~ f . The vertex correction leads to further reduction of the effective charge. Finite contributions 
to dielectric screening are expected in all orders of perturbation theory. 
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I. INTRODUCTION AND MOTIVATION 

Graphene is a two-dimensional (2D) allotrope of car- 
bon on a honeycomb lattice with one electron per 7r or- 
bital (half- filling) . Its bare electronic spectrum is de- 
scribed in terms of a linearly dispersing, massless, chi- 
ral Dirac field (\& p ). Since its isolation a few years 
agoi it was realized that graphene displays a number 
of unique properties that are at odds with the stan- 
dard theory of metals* 1 -^ One of the most important 
unresolved questions in graphene is the role of electron- 
electron interactions 4£&L£ & 10 i 11 i 12 Even though, due to 
the vanishing of the density of states at the Fermi energy, 
the electron-electron interactions are expected to remain 
unscreened and strong, it is not clear what is their influ- 
ence in the properties of graphene. 

In the present work we study the influence of the 
electron-electron interactions on the static dielectric 
function of graphene at half filling. We perform cal- 
culations to one order beyond the conventional random 
phase approximation (RPA) vacuum polarization bub- 
ble, by including self-energy and vertex corrections into 
the polarization loop. Our main finding is that vertex 
contributions become important in the coupling regime 
a ~ 1, which in turn means that other non-RPA contri- 
butions should also be included. We were mainly moti- 
vated by the question whether the interactions can sig- 
nificantly affect the screening properties. This issue is 
particularly relevant in graphene for two reasons: (1.) 
The effective coupling constant a (see the precise defi- 
nition below) in graphene is large a ~ 1, and thus in- 
teractions are expected to be generically important, and 
(2.) Despite of the above, to the best of our knowledge, 
no clear signatures of interaction effects have been ob- 
served so far in graphene. For example measurements of 
the compressibility^ have not detected electron correla- 
tion effects. In addition, screening of external charged 
impurities introduced in graphene is also expected to 
be sensitive to interaction effects, at least on theoretical 
level) 14 : 15 ' 16 : 17 and could be relevant for interpretation of 
recent experiments on charged impurity scattering^ It 
is thus generally important to investigate the problem 
of how the correlations affect the effective charge of the 
carriers in graphene, which is determined by the vacuum 



polarization. We will assume that graphene at half filling 
(i.e. when the chemical potential crosses the Dirac point) 
remains a homogeneous gas of quasiparticles, which in it- 
self is not necessarily an innocent assumption due to the 
possibility of puddles, ripples, etc4 However we assume 
that the system is homogeneous as the importance of the 
above effects is still unsettled. 

Our starting point is the low-energy Hamiltonian of 
graphene which can be written as (we use units such that 

h = i), 

H = H*p( u IpI^p- »vo)V P + Hi, (1) 
p 

where v is the Fermi velocity, fj, is the chemical potential 
away from half- filling, &o = I is the 2x2 identity matrix, 
<7p = a ■ p/|p| = (a x p x + ayPy^p^ 1 , and a x ,a y are 
Pauli matrices. The first term in the Hamiltonian ([T]) 
reflects the effective Lorentz invariance that exists in the 
non-interacting problem at low energies and gives rise to 
bizarre electronic behavior analogous to the one found in 
quantum electrodynamics (QED). 19 

In Eq.JT]), Hj is the electron-electron interaction, 

ff/=2&pV^P' ™p =J2K+p*<i> ( 2 ) 
p q 

where 



is the Fourier transform of the Coulomb potential in 2D. 
The relative strength of the Coulomb interactions to the 
kinetic energy is determined by graphene's "fine struc- 
ture constant" a = e 2 /v. Unlike QED, the Dirac fermion 
velocity is much smaller than the speed of light, c, and 
hence the Coulomb field can be treated as instantaneous 
(v w 10 6 m/s). As a result, the Coulomb interaction 
breaks the Lorentz invariance of the problem leading to 
fundamental differences between the graphene problem 
and QED. From now on we absorb the dielectric con- 
stant of the medium s into the definition of the effective 
charge e. For example in the typical case of graphene 
on a Si02 substrate with dielectric constant e ps 4, we 
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have e 2 = 2e^/ (1 + e), where eo is the charge of the elec- 
tron. Keeping in mind that e^/v w 2.2, one then finds 
the coupling constant a « 0.9. — Nevertheless, even in 
this situation when the relation a < 1 is not strictly 
satisfied, perturbation theory is expected to give a good 
indication for the behavior of physical quantities. 

The rest of the paper is organized as follows. Section 
II deals with corrections to the polarization loop arising 
from the dressing of the electron propagators. In the 
Section III the true interaction (correlation) insertion, 
the vertex correction, is examined. Section IV contains 
our conclusions. 



II. SELF-ENERGY CORRECTIONS TO THE 
POLARIZATION 

We concentrate on the most interesting case of zero 
Fermi energy (/i — > 0), when the low-energy physics is 
controlled by the proximity to the Dirac point. 

The free Dirac Green's function is 

G(k, w) = I - n + • ( v W 

woo — v|k|ak + «croO + sign(ijj) 

The interaction effects lead to the dressed Green's func- 
tion G _1 — > G -1 — £, where the self-energy £ is a sum of 
two terms with different matrix structure: E = Eo + S, 
Eo oc <To,E oc <Tk- At Hartree-Fock (HF) level (first 
order in a) a divergent contribution appears, due to 
the long-range nature of the Coulomb interaction^ where 
A^l/a^fcisan ultraviolet cutoff (a is the lattice 
spacing). One finds 



£«(|k|) = (a/4) Hk|)<7 k ln(A/|k|), 



(5) 



which implies that the effective velocity changes v — > 
v(l + (a/4) ln(A/|k|)), and grows without bound at low 
energies |k|/A — * 0. This should in principle lead to 
anomalies in thermodynamic and spectral properties of 
graphene4^i From theoretical viewpoint, most impor- 
tantly, the single logarithmic behavior was found to per- 
sist to second order of perturbation theory as wellj& and 
consequently this is expected to be the case to all orders, 
reflecting the fairly simple (at least at weak-coupling) 
renormalization structure of the theory. 

We now turn to the calculation of the static polar- 
ization, II(q) = II(q, u) = 0). The frequency variable in 
II(q) is omitted from now on. The bare polarization bub- 
ble (without any interaction lines in the loop) is known 
to be 



n<°>(q) = "*£/ ^Tr{G(k,u,)G(k + q,u,)} 



= -|q|/(4«). 



(6) 



From now on the trace stands for summation over spin 
(s), valley (v) and pseudospin (Pauli matrix a) indices, 
i.e. 



Tr = ^Tr ff = 4Tr 



(7) 
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FIG. 1: First order interaction corrections to the polarization 
bubble: (a.) self-energy correction, (b.) vertex correction. 
The explicit form of these diagrams is given by equations (9) 
and (19). The wavy line represents the Coulomb interaction, 
Eq. (3). 



The momentum sums are performed as = 
Jd 2 k/(2n)*. 

Next, we calculate the bubble dressing due to the 
electron-electron interactions to first order in a. The 
two diagrams at this order are shown in Fig. [T] We write 
the total polarization as 



n( q ) = n(°) ( q ) + n« ( q ) + n« ( q ) 



(8) 



where IlW and II^ 2 ) stand for the contributions of 
Fig. QJa) and Fig. QJb), respectively. 

The self-energy dressing of Fig. [Ha) can be written as 
(the factor of 2 originates from the two possible inser- 
tions) 

nW(q) =-2iJ2 I ^Tr{G(k,c)G(k + q,c)G(k, W ) 



xE(k, W )}. 



(9) 



At lowest order, the self-energy is simply the Hartree- 
Fock one, meaning that in ^ we replace 



E(k,w) -► E (1) (k) 



X / ^G(p, Wl )V- k _ p . (10) 



The large logarithm present in EW(k) at low momenta, 
Eq. ([5]), is expected to appear also in some form in 

n«(q). 

Let us define the following quantities which appear in 
our results from now on 

A(k,q)^l- k ; i (k + q) , fc = k/|k|, (11) 



£7(k,q) = t;(|k| + |k + q|; 



(12) 



After performing the frequency, and then momentum in- 
tegrations in © , and using the self-energy from (IL0|) , we 
obtain 



nW(q) =4^Vk-p(k-p) 

k,p 



A(k,q) 
[S(k,q)] 2 



a_ |q| 
16 v 



ln(A/|q|), A/|q|»l. (13) 
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FIG. 2: First order interaction correction to the vertex T(k;q), 
given by Eq. (|T¥|) . 

This result means that the large logarithm in I]W(q) 
simply reflects the renormalization of the Fermi veloc- 
ity, i.e. this divergence is not independent, but can be 
simply reabsorbed into the velocity by replacing v — > 
u(l + (a/4)ln(A/|q|)) in the one-loop result E[W>(q) = 
— |q|/(4u). Due to the simple logarithmic structure of 
the theory this is expected to hold to all orders of per- 
turbation theory, i.e. all self-energy corrections lead to a 
replacement of the coupling a in all final expressions with 
the "running" coupling a(q), accounting for the velocity 
renormalization. We therefore assume that the velocity 
renormalization procedure is performed in all higher or- 
der diagrams. At finite small chemical potential, /i <C A, 
which is the case in any realistic experimental situation, 
the divergence is cut-off, ln(A/|q|) — * ln(A//i). Due to 
the slow variation of the logarithmic function or possi- 
ble other factors (such as strong dielectric screening), no 
significant variation of the velocity has been found in 
experiment. 3 

An interesting effect, related to the interaction con- 
tribution nW(q), Eq. (|13[) . was recently discussed in 
Ref. [l6| within the RG approach. Our calculation, lead- 
ing to Eq. (|13|) . provides an explicit perturbative con- 
firmation of the RG results. If we imagine an external 
Coulomb impurity with charge Z > 0, probing the po- 
larization of the vacuum, then the induced charge den- 
sity, in momentum space, is Pi nc t(<l) = ZV c JI(q L ). Here 
n(q) = n(°)(q) + nW(q). While the first term leads 

to induced charge ^^(q) = —Za(ic/2), localized in real 
space at the impurity site and with a screening sign, the 
interaction term pi„j(q) = Za 2 (n/8) ln(A/|q|) has an 
opposite sign and decays as an inverse power law (1/r 2 , 
with logarithmic corrections). This peculiar behavior 
simply reflects, however, the renormalization (increase) 
of the Fermi velocity at low momenta, which leads to 
suppression of screening at large distances. 



III. VERTEX CORRECTIONS TO THE 
POLARIZATION BUBBLE 

We proceed with the calculation of the vertex correc- 
tion in Fig. HKb). Before evaluating this expression, it 



is useful to examine the (possible) singularity structure 
separately in the vertex operator, shown in Fig. [21 For 
simplicity we use the notation T(k;q) which stands for 
the more conventionally used form F(fc, k + q;q), where 
the variables denote both frequency and momenta (q is 
the bosonic momentum/frequency). 

Since the Coulomb interaction is non-retarded, we have 
the simple expression 

^G(p, Wi)G(p + q, wi + u)V^ p . 

(14) 

After evaluating the frequency integral, the result is a 
sum of an off-diagonal and diagonal parts (with respect 
to the Pauli matrix indexes), f = + f ( d \ More 
explicitly, 

f<«>(k;q,*) ccX> k _ p (4„ - » p+q ) ( E{p \ )+ul 

-HS5F=)- ,15 » 

We are interested only in the zero frequency limit (and 
only in the real part of T, since the imaginary part does 
not contribute to the polarization). In this case the off- 
diagonal piece vanishes identically, r(°'(k;q, uj = 0) = 0. 
This is expected to be the case since the Coulomb inter- 
action is diagonal and thus the vertex can not generate 
a static contribution with a different matrix structure. 
On the other hand the diagonal part is finite in the same 
limit 

f ^(k; q,^0)«^ Vk- P (l - a p( 7 p+q )-^i-^. (16) 

An explicit evaluation shows that r( d '(k;q, u = 0) does 
not have any divergent contributions. For example at 
|k| « |q|, 

f (d) (k;q,0) oc a const., |k| < |q|, (17) 

while in the opposite limit 

fW(k;q,0)<xa(|q|/|k|), |q| « |k| (18) 

For our purposes the exact formulas are not important 
(we also do not show the dependence on the angle be- 
tween k, q) ; our main conclusion at this stage is that the 
vertex does not have any divergent parts. We have also 
examined diagrams of higher order, such as "ladder" and 
"crossed" ladder vertex corrections, and have found that 
all of them are finite. Therefore the vertex insertions 
into the polarization function are expected to give a fi- 
nite contribution to that quantity, and below we evaluate 
the lowest order vertex correction numerically. 

It is clear that a Ward identity relating divergent con- 
tributions in the self-energy and in the vertex does not 
hold here, unlike conventional QED where Lorentz (and 
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gauge) invariance guarantees cancellation between ver- 
tex and self-energy corrections^ and charge is renormal- 
ized only through simple polarization loops in the photon 
propagator. On the other hand in graphene, where the 
only non-trivially renormalized quantity is the velocity v, 
both the polarization operator and the vertex operator 
do not show any independent divergencies. 
The diagram of Fig. [T]Jb) now reads 

n (2) (q) = I ?Tr{G(k, W )f (k; q, 0)G(k + q, 

k J n 

(19) 

where the full expression for T from Eq. (jT4j) should be 
used. An explicit calculation, starting by evaluation of 
the energy integrations, leads to the result 



k, P 

Taking into account: 

1 

CkCp 



£(k,q)£(p,q) 



(20) 



W|p| 



(k.p + w 3 (kxp) z ), (21) 



where (p x q) stands for a vector product, we then arrive 
at the final formula 



lK 2 )(q) =-2£ 

k,p 



Vk-p 
£?(k,q)£(p,q) 



{A(k,q)A(p,q) 



(P x q) 2 (k x q) z 
|p||k||p + q||k + q| 



(22) 



It is clear on dimensional grounds that n( 2 )(q) is linear 
in |q|. This is in fact the case for polarization diagrams in 
all orders of perturbation theory. The four-dimensional 
integrals, appearing in (|22|) , cannot be evaluated analyti- 
cally. We have found, as expected in light of our previous 
discussion of the vertex function, that the expressions 
converge in the ultraviolet limit. After computing the 
integrals numerically, we obtain the following result for 
the combination VqII^ 2 ^(q), which appears in the dielec- 
tric function, 



2ne 2 



n (2) (q) = -0.53a 2 . 



(23) 



Adding also the one- loop RPA result, we have finally 
(where £ is the static dielectric constant, defined by the 
formula below) 



V r 



off. 



V a 



q 1 - y q n( q ) £ ' 



£ = 1 



0.53 a 2 + 0(a 3 



(24) 



(25) 



We conclude that, at a ~ 1, the vertex correction 
is more than 30% of the one-loop result. It also has a 



screening sign, i.e. it reduces the effective charge. One 
also expects that finite contributions will appear to all or- 
ders in a. However, resummation of perturbation theory 
by simple means seems impossible, as the contributions 
in question are finite and accumulate over a wide range 
of momenta in the corresponding diagrams (rather than 
within a specific integration window, from where diver- 
gent parts typically originate, and thus can be easily col- 
lected). Even though the vertex contribution is a sizable 
one, two remarks are in order: (1.) It does not change 
drastically the structure of the theory, apart from con- 
tributing towards further screening of the interactions. 
(2.) The fact that perturbation theory is used with the 
intention of being applied at a rather strong coupling is 
in itself questionable. Nevertheless, perturbation theory 
provides a clear indication that a significant contribu- 
tion to screening exists beyond the conventional one-loop 
RPA result. On the other hand in the weak-coupling 
regime, a < 1, RPA is parametrically well justified as 
far as the static polarization properties are concerned 
(although the RPA is not justified for the self-energy.—) 



IV. DISCUSSION AND CONCLUSIONS 

It is also useful to compare our results to the situation 
in ordinary metals with a finite Fermi surface. Certain 
approximations are typically used to account for vertex 
corrections, such as the Hubbard form of the dielectric 
function. When extrapolated to low momentum, the ver- 
tex contribution tends to decrease the screening length, 23 
i.e. it reduces further the range of the interactions. Nat- 
urally in graphene, where the screening length is infinite 
(for the case of zero chemical potential considered here), 
the vertex correction affects directly the effective charge, 
without changing the shape of the Coulomb potential. 

Finally we mention two recent related works, dis- 
cussing interaction effects, that appeared while the 
present manuscript was being prepared. In Ref. I24L the 
effect of self-energy and vertex corrections to lowest or- 
der (a) on the minimal conductivity in graphene was dis- 
cussed, with the conclusion that the corrections is of or- 
der 1%. Dynamical polarization properties were studied 
in Ref. l25l . where the vertex diagrams were found to have 
logarithmically singular contributions near the threshold 
lu ~ vq, leading to the possibility of a plasmon mode. 

In summary, we have shown that vertex corrections 
can have sizable effect in the static vacuum polariza- 
tion diagrams in the regime of strong coupling, while for 
small coupling their importance diminishes parametri- 
cally. The self-energy corrections are naturally absorbed 
into the renormalization of the Fermi velocity. The non- 
RPA vertex diagram at lowest order of perturbation the- 
ory was found to decrease the effective charge, meaning 
that in principle correlation effects at higher order must 
also be taken into account. Thus the ultimate asymptotic 
behavior of the static polarization function for a ~ 1 re- 
mains an open problem. 
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